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Abstract 

Background: Malaria control programmes across Africa and beyond are facing increasing insecticide resistance in the 
major anopheline vectors. In order to preserve or prolong the effectiveness of the main malaria vector interventions, 
up-to-date and easily accessible insecticide resistance data that are interpretable at operationally-relevant scales are 
critical. Herein we introduce and demonstrate the usefulness of an online mapping tool, IR Mapper. 

Methods: A systematic search of published, peer-reviewed literature was performed and Anopheles insecticide 
susceptibility and resistance mechanisms data were extracted and added to a database after a two-level verification 
process. IR Mapper (www.irmapper.com) was developed using the ArcGIS for JavaScript Application Programming 
Interface and ArcGIS Online platform for exploration and projection of these data. 

Results: Literature searches yielded a total of 4,084 susceptibility data points for 1,505 populations, and 2,097 resistance 
mechanisms data points for 1,000 populations of Anopheles spp. tested via recommended WHO methods from 54 
countries between 1954 and 2012. For the Afrotropical region, data were most abundant for populations of An. 
gambioe, and pyrethroids and DDT were more often used in susceptibility assays (51.1 and 26.8% of all reports, 
respectively) than carbamates and organophosphates. Between 2001 and 2012, there was a clear increase in 
prevalence and distribution of confirmed resistance of An. gambioe s.l. to pyrethroids (from 41 to 87% of the mosquito 
populations tested) and DDT (from 64 to 91%) throughout the Afrotropical region. Metabolic resistance mechanisms 
were detected in western and eastern African populations and the two kdr mutations (L1014S and L1014F) were 
widespread, for An. funestus s.l., relatively few populations were tested, although in 2010-2012 resistance was reported 
in 50% of 10 populations tested. Maps are provided to illustrate the use of IR Mapper and the distribution of insecticide 
resistance in malaria vectors in Africa. 

Conclusions: The increasing pyrethroid and DDT resistance in Anopheles in the Afrotropical region is alarming. Urgent 
attention should be afforded to testing An. funestus populations especially for metabolic resistance mechanisms. IR 
Mapper is a useful tool for investigating temporal and spatial trends in Anopheles resistance to support the pragmatic 
use of insecticidal interventions. 
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Background 

Malaria remains one of the major disease burdens globally 
with over 207 million cases and 627,000 deaths estimated 
in 2012, mainly in children under 5 years old and predo- 
minantly in Africa [1]. The Democratic Republic of Congo 
(DRC) and Nigeria account for 40% of the total estimated 
global deaths, and combined with India account for 40% 
of the total estimated cases. Malaria is considered to 
be one of the major contributors to poverty and the 
estimated annual cost to economies of the African conti- 
nent ranges from under 0.5% to almost 9% of GDP [2]. At 
1% of the total African GDP of US$ 1,184 trillion, this 
translates into US$ 12 billion per annum. In addition to 
the scale of mortality and the loss of productivity due to 
illness, malaria also has devastating effects on cognitive 
development in children surviving the disease, leaving 
many disabled for life [3,4]. 

The dominant mosquito species responsible for the 
transmission of malaria parasites in Africa are mainly 
Anopheles gambiae s,s,, An, coluzzii (formerly gambiae 
s.s. M form), An. arabiensis and An, funestus s.s. that are 
widespread over tropical and subtropical Africa, although 
An. arabiensis prefers drier habitats and An, coluzzii is re- 
stricted to west-central Africa [5,6]. The adult behaviour 
and larval biology of these species are different, which im- 
pacts the ease with which each can be controlled. Adult 
An, gambiae s,s„ An, coluzzii and An, funestus s,s, prefer 
feeding on humans and resting inside human habitations, 
while An, arabiensis will feed on either humans or cattle 
and rest indoors or outdoors making this vector more dif- 
ficult to control [5-7]. Anopheles funestus prefers to breed 
in swampy, well-vegetated, permanent water bodies, while 
An, gambiae s,s„ An, coluzzii and An, arabiensis can be 
found in temporary rain pools, hoof prints around the 
edges of dams and pans, and also in rice paddies [5-7]. 
Larviciding and/or larval habitat modification is therefore 
not always simple or feasible. 

The most widespread and effective approach to control- 
ling malaria vector mosquitoes is through the use of in- 
secticide treated bed nets (LLINs) and indoor house 
spraying with residual insecticides (IRS). By far the most 
common insecticides used for house spraying are pyre- 
throids and this class of insecticide is the only one cur- 
rently recommended by the World Health Organization 
(WHO) for treatment of bed nets [1]. It is therefore not 
surprising that malaria vectors have developed resistance 
to pyrethroids throughout the African continent [8-10]. 
Three other classes of WHO-recommended adulticides 
(organophosphates, carbamates and organochlorines) have 
also been used in IRS to differing extents throughout 
Africa, and resistance has been detected in Anopheles spp. 
to all three [11]. Insecticide resistance data are needed to 
inform vector control policy, particularly in the context of 
the Global Plan for Insecticide Resistance Management 



(GPIRM) that seeks to preserve or prolong the effective- 
ness of vector control interventions [11]. 

For optimal utility in informing policy, vector data need 
to be accurate, up-to-date, easily accessible and inter- 
pretable at operationally relevant scales. A number of pa- 
pers have produced maps of malaria vector species on a 
country level [12-15] and across the African continent 
[6,9,10,16-18], mostly dealing with the distribution of the 
vectors but very few with insecticide resistance [8-10]. 
Methods used range in complexity, from simple presence/ 
absence plotting on a map [8,9,12,16,17] to more sophisti- 
cated predictive models [17,19-23]. The difficulties facing 
interpretation of historical species occurrence records 
such as poor data coverage and taxonomic ambiguity of 
species also affect the mapping of insecticide resistance. 
Other issues facing the consolidation of insecticide resis- 
tance data include a bias towards reporting of susceptibil- 
ity results only when resistance is subsequently detected 
[24,25], insufficient information to determine if standard 
protocols were adhered to, incorrect or missing geo- 
coordinates for collection sites, and a proliferation of pub- 
lished data for sentinel sites favoured by research groups 
while large volumes of data for other regions remain 
unpublished [9]. Disparities in criteria for reporting resis- 
tance also complicate matters - this is of particular 
relevance given a recent increase in the threshold for 
reporting resistance from 80% to 90% mortality in WHO 
susceptibility tests [26] . 

Static maps in publications quickly become outdated due 
to rapid changes in insecticide resistance status being re- 
ported in the current relatively high volume of publications 
on the topic. Until recently there was no consolidation of 
all historical and up-to-date information on insecticide re- 
sistance in malaria vectors. While there have been a num- 
ber of initiatives to create online repositories specifically for 
such data, including within the Mapping Malaria Risk in 
Africa (MARA) [9] and VectorBase [27] databases, these 
platforms did not include quality control for uploaded data 
or failed to present summaries of recent data in a format 
appropriate for informing vector control decisions. 

In this paper we introduce the free online geospatial ap- 
plication IR Mapper (www.irmapper.com), which facilitates 
the exploration and projection of worldwide insecticide 
resistance data. We also provide an update on the status of 
insecticide resistance in the main malaria vectors in Africa 
including an historical review based on published reports 
from insecticide susceptibility tests and mechanisms inves- 
tigations. Data are reviewed with respect to updated WHO 
guidelines refining the definition of "resistance" [26]. 

Methods 

Insecticide resistance database 

A systematic search of the published, peer-reviewed litera- 
ture using online scientific bibliographic databases was 
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performed using key words in English and French. Data- 
bases included PubMed [28], Web of Science [29], ar- 
chives of Malaria World [30], Google Scholar [31], and the 
Armed Forces Pest Management Board [32]. Similarly, ar- 
chives of 10 top journals in the field of Anopheles and in- 
secticide research (Malaria Journal, Parasites and Vectors, 
PLOS ONE, Medical and Veterinary Entomology, Journal 
of Medical Entomology, Tropical Medicine and Inter- 
national Health, American Journal of Tropical Medicine 
and Hygiene, Transactions of the Royal Society of Tropical 
Medicine and Hygiene, Journal of Vector Ecology, Journal 
of Vector Borne Diseases) were searched. Key words used 
in singular or combination and in English and French 
were acetylcholinesterase. Anopheles, carbamate, DDT, 
esterase, GST, insecticide, kdr, pyrethroid, metabolic, 
mono-oxygenase, organochlorine, organophosphate, oxi- 
dase, resistance, resistant, target site, and names for indi- 
vidual countries. The search exercise was completed for 
all 54 countries of the African region. This process was 
similarly conducted for other countries outside Africa 
(17 to date). 

Reference sections of all relevant located articles were 
also reviewed to identify additional literature. The data- 
set was augmented with three extra unpublished sources 
of information. The first was the African Network for 
Vector Resistance (ANVR) data extracted from IRBase 
[33]. The second was a summary report by the Presi- 
dents Malaria Initiative (PMI) including WHO suscepti- 
bility test and kdr data from 18 countries [34]. The third 
source was unpublished reports of WHO susceptibility 
tests conducted using standard procedures on Anopheles 
populations from Cote d'lvoire (2008, 2010), Democratic 
Republic of Congo (2007, 2011), Mali (2000, 2003, 2005, 
2006, 2007, 2011) and Tanzania (2004, 2012) (Richard 
Hunt, unpublished data). 

Data were extracted into Microsoft Excel 2010 (Redmond, 
USA) datasheets. For each insecticide susceptibility or 
mechanisms test conducted, the following were recorded: 
country name, locality name, GPS coordinates (latitude, 
longitude) of locality, mosquito collection period, vector 
species or species complex/group tested. Where data on 
mosquito collection period were not reported in the publi- 
cation, efforts were made to contact the authors to obtain 
this information. In the event that no response was re- 
ceived from the authors, the year of publication was taken 
as the year of mosquito collection. For WHO suscepti- 
bility tests on adult mosquitoes, the following were 
recorded: insecticide type, insecticide class. Insecticide 
Resistance Action Committees (IRAC) mode of action 
classification and code, insecticide dosage, number of 
mosquitoes exposed per assay, measured percentage mor- 
tality corrected for controls, and susceptibility status. 
Status was assigned based on recently-revised WHO cri- 
teria [26], where: mortality <90% = confirmed resistance; 



90-97% = possible resistance (with presence of resistant 
genes to be confirmed); 98-100% = susceptible. Where 
mortality rates were reported in range format, the average 
of the highest and lowest values was used to assign resis- 
tance status. 

For molecular or biochemical mechanisms data, the fol- 
lowing were recorded: test method used, mechanism 
tested for, number of mosquitoes used per assay, outcome 
of the assays (detected/ not detected), frequency of muta- 
tions. No assumptions were made in the data abstraction, 
with all reported data accurately reflecting the level of 
detail given in the data source. Geo-referencing was con- 
ducted in decimal degrees format using the set of coordi- 
nates provided in publications. Where coordinates were 
not listed, these were determined by locality names via 
Geonames [35] and failing that, Google Earth [36]. Where 
the administrative unit name was given only, this was in- 
cluded in the database with coordinates assigned based on 
the centre point of the unit; administrative units ranged 
from districts or counties to sub-districts or sub-counties 
depending on the particular publication. Data from local- 
ities for which coordinates could not be located were 
excluded from the database. 

Data extractions primarily focused on publications with 
the assumption that published data were quality assured. 
However, some of the extracted data did not adhere to 
standard WHO protocol. For instance, in WHO suscepti- 
bility tests, assays were conducted using insecticide papers 
with non-standard concentrations of insecticide, or fewer 
than the recommended minimum of 100 mosquitoes were 
tested to derive a mortality rate. For this reason, such 
major parameters influencing study outcomes were also in- 
cluded in the database for presentation on the IR Mapper 
interface and were considered during the review process, 
as outlined below. 

Foflowing extraction, data were subjected to a two- 
level checking procedure, with the first level check 
conducted by a different abstractor to ensure an inde- 
pendent assessment of the assembled data. Data were 
cross-checked with any repeat records removed, espe- 
cially for instances where unpublished PMI or ANVR 
data had subsequently been published. AU aspects of the 
data were reviewed to ensure the information had been 
correctly collated and the sites geo-referenced accur- 
ately. The second level check incorporated suggested 
changes with an emphasis on geo-referencing. Third- 
level checks were implemented using the interface out- 
lined below to identify: (i) inconsistent or non-standard 
spellings within fields that would affect filter query sum- 
maries; (ii) blanks in mandatory data fields; (iii) any 
geo-locations that fell outside the specified country 
boundaries or in the sea or other major water body; and 
(iv) any localities with similar names but different coor- 
dinates in order to harmonize naming/location data. 
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Literature searches and data extraction for the purpose 
of this review were concluded on 15 May 2013. How- 
ever, the intention is for the process described above to 
continue on a monthly basis to ensure regular update 
of the database accessible via the IR Mapper online 
interface. 

The database contained information extracted from 224 
individual publications from 60 peer-reviewed journals. 
There were relatively few publications up until the 
1990s, with a rapid increase in both insecticide suscepti- 
bility and resistance mechanisms reported after 2005 
(Figure 1). For the first four months of 2013, there were 
14 publications that included data on susceptibility and 
16 on mechanisms (18 in total). It is likely, therefore, 
that the number of publications for 2013 may exceed 
that of previous years (maximum of 26 in 2012). The 
database includes a total of 4,084 insecticide suscepti- 
bility tests for 1,505 Anopheles populations collected 
from 1954 to 2012, and 2,097 resistance mechanisms 
tests for 1,000 populations collected from 1987 to 2012. 



mechanism/s and dates of mosquito collections. A time 
filter enables users to view data by single or multiple year 
increments based on the start date of mosquito collections. 
Plasmodium falciparum and Plasmodium vivax ende- 
micity maps [37,38] provided by the Malaria Atlas Project 
[39] are incorporated as optional layers. User-generated 
map images can be printed or saved. 

Extended functionality allows users to download a stan- 
dardised Excel template, add their own geo -referenced 
data and temporarily visualise these data on the mapping 
interface. Filters can then be applied to users' own and 
existing data, and if desired the latter can be hidden via 
the left hand data layer menu. Therefore, this provides a 
facility for users to visualise (and print) all or part of their 
own unpublished data set either independently or along- 
side the published data set. Users' data are not added to 
the online database and are cleared once the browser is 
closed. This ensures that all data available on the IR 
Mapper interface have undergone the quality control and 
assurance processes outlined above. 



Online mapping interface 

An online interactive mapping platform, IR Mapper (www. 
irmapper.com), was developed to allow users to project 
geo-referenced data residing in the database (Figure 2). 
The web application was built on the ArcGIS Application 
Programming Interface for JavaScript. Outcome data from 
WHO susceptibility tests (confirmed resistance, possible 
resistance, susceptible) and mechanisms assays (detected, 
not detected) are displayed. Details of assays (including in- 
secticide dosages and numbers of mosquitoes tested) and 
links to data sources are provided in pop-up boxes asso- 
ciated with individual point data. The mapping application 
facilitates filtering of data via user-specified criteria for 
Anopheles species, insecticide class/es or type/s, resistance 



Review of insecticide resistance status for the Afrotropical 
Region 

A sub-set of the full database described above and avail- 
able on IR Mapper was extracted for the purpose of this 
review. Date of acceptance for publication was used as 
the cut-off, with data from articles accepted after 30th 
April 2013 excluded. PMI [34] and ANVR [33] unpub- 
lished datasets were excluded on the basis that data 
ownership resides with the institute/s responsible for 
collecting and reporting the data for their region. Results 
from WHO susceptibility tests with insecticidal papers 
that were not of standard discriminating dosages accor- 
ding to WHO [26] were excluded from the review, with 
the exception of permethrin and deltamethrin for which 
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Figure 1 Number of publications from which data were extracted on insecticide susceptibility and resistance mechanisms, up to 30th 
April 2013. 
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Figure 2 IR Mapper online user interface showing Anopheles spp. pyrethroid resistance (red dots) and detected elevated mono- 
oxygenase and esterase activity plus kdr mutations (black dots) from 2001 - 2013 (as at 06/12/2013) along with Plasmodium falciparum 
endemicity estimates for 2010. 



different standard dosages (0.25% and 0.025%, respec- 
tively) applied from 1981 [40] until 1998 [41] (Figure 3). 
Since discriminating dosages have yet to be defined for 
alpha-cypermethrin and chlorpyrifos-methyl, suscepti- 
bility data for these insecticides M^ere excluded from the 
data set. 

Scope was limited to tests on adults of the main African 
malaria vector species of the An. gambiae complex and 
An, funestus group, with reviews conducted based on four 
groupings: a) all species of the An, gambiae complex; 
b) An, gambiae s,s,; c) An, arabiensis; and d) all species of 
the An, funestus group. For the purpose of the review. 
An, gambiae s,s, (previously the S molecular form) and 
An, coluzzii (formerly An, gambiae s,s, M form) were 
grouped together since the nomenclature change is recent 
[5] and thus is yet to be consistently adhered to in publica- 
tions. A population was defined as any species or species 
complex/group collected from one reported location and 
at one reported time period. The duration of collections 
varied by study, and was not always specified; these may 
have varied from days to months. For review purposes, ag- 
gregation of data by year was based on the start year of 



the collection period. Predicted distribution maps by vec- 
tor species were sourced in GIFF format [42]; individual 
species maps were used for An, gambiae s,s, and An, 
arabiensis or maps of individual species were merged to 
create maps by An, gambiae complex or An, funestus 
group. All maps for the purposes of the review were 
created in ArcGIS for Desktop 10.1 (Esri, Redlands, CA). 

Results 

The example screenshot of the IR Mapper online user 
interface shown in Figure 2 depicts recorded pyrethroid 
resistance (red dots) and detected elevated mono- 
oxygenase activity, elevated esterase activity, and/or kdr 
mutations (black dots) for Anopheles spp. populations 
collected from 2001 to 2013 (as at 06/12/2013). These 
point data are overlaid with a map layer showing esti- 
mated levels of Plasmodium falciparum malaria endemi- 
city within the limits of stable transmission for 2010 [37]. 

Additional maps show the capacity for data residing in 
IR Mapper to be used to illustrate more detailed aspects 
such as temporal and spatial distribution in insecticide 
resistance and the known mechanisms underpinning 
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Figure 3 Publication dates of Icey protocols and procedures used in insecticide resistance evaluations for Anopheles. Ace-IR: insensitive 
acetylcliolinesterase Gl 19S mutation; CDC: Centres for Disease Control and Prevention; GSTs: elevated glutathione s-transferase; kdr. L1014S and L1014F 
target site mutations; MR4: Malaria Research and Reference Reagent Resource Centre; PCR: polymerase chain reaction; PCR-RFLP: PCR restriction 
fragment length polymorphism; rdl: dieldrin resistance-conferring chloride ion channel mutation; RT-PCR: reverse transcriptase PCR; SNP-PCR: single 
nucleotide polymorphism PCR; WHO: World Health Organization. 



resistance for each of the major malaria vectors or 
species groups. 

Susceptibility bioassays 

Insecticide susceptibility data were recorded for 28 coun- 
tries and resistance mechanisms data for 31 countries out 
of the 54 in the African region. Data on susceptibility were 
most numerous for Benin, Cameroon and Mozambique 
(92, 81 and 79 populations, respectively) and for mecha- 
nisms were the most numerous for Burkina Faso, Benin 
and Kenya (166, 128 and 102 populations, respectively). 
Both data types were most often reported for Anopheles 
gambiae si (716 populations) with significantly fewer re- 
ports for members of the An, funestus group (123 popu- 
lations); bioassay data were rarely reported for sibling 
species indicating how seldom molecular species differen- 
tiation was conducted for bioassay specimens (Additional 
files 1 and 2). The number of populations tested seemed 
to decrease in 2010-2012, but this is likely because results 
obtained for collections from this period are yet to be pub- 
lished. Across species and time, pyrethroids were generally 
the most tested insecticide class (51.1% of all test reports), 
followed by organochlorines (26.8%), carbamates (12.9%), 
and organophosphates (9.2%). Testing for kdr mutations 
was most common for An, gambiae s,l, with testing for 
metabolic mechanisms seldom conducted. For the period 
2001-2012, the majority of insecticide susceptibility (78%) 
and resistance mechanism (73%) tests were reported for 
Anopheles populations from the 20 top malaria burden 
countries in Africa although relatively few data were avail- 
able for some key countries such as DRC (Additional file 3). 



During 2001-2012, resistance was detected to at least 
one insecticide in 78.5% (472/601) oiAn, gambiae s,l, pop- 
ulations from 27 of the 28 countries tested (with suscepti- 
bility only detected in Guinea Bissau) and in 43.6% (34/78) 
of An, funestus from 11 out of 13 countries tested (with 
susceptibility only detected in Burundi and Tanzania). 
Resistance to two classes of insecticides was detected in 
219 populations oiAn, gambiae s,l, and 11 oi An, funestus 
while resistance to three classes was found in 40 and 2 
populations of these species, respectively. Resistance to in- 
secticides from all four classes was detected in individual 
An, gambiae s,l, populations from five localities in Burkina 
Faso, three in Cote d'lvoire and one each in Mali and 
Sudan. The highest numbers of reports of confirmed re- 
sistance since 2001 were for Anopheles populations from 
Benin (151), Cameroon (104) and Nigeria (94). 

The prevalence of An, gambiae s,l, resistance to pyre- 
throids and DDT increased between the periods 
2001-2003 and 2010-2012: for pyrethroids from 41% to 
87% and for DDT from 64 to 91% (Figure 4, Additional file 
4). In An, funestus y pyrethroid resistance prevalence in- 
creased from 26% of 31 populations pre-2001 to 50% in 
2010-2012 for 10 populations; a significant drop in 2004- 
2006 was observed though it is important to note that 
populations tested in the previous time period (2001- 
2003) were all from a single country (Mozambique). The 
prevalence of confirmed resistance of An, funestus to 
DDT was between 0% and 28%, though the lowest figure 
was for only two populations (Figure 4, Additional file 4). 
Similar examinations indicated no dramatic increases 
in the prevalence of carbamate or organophosphate 
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Figure 4 Proportion of populations of An. gambiae s.l. and An. funestus s.l. with confirmed resistance using standard WHO susceptibility 
tests for pyrethroids and organochlorines. Note that populations are accounted for more than once when tested with both insecticide 
classes. Italics indicate the number of populations tested. 



resistance in either species group. There have been no re- 
ports of resistance of An. funestus to organophosphates 
since 1978. 

Resistance mechanisms 

Investigations of resistance mechanisms in An. gambiae 
have largely focussed on target site mutations {kdr and 
Ace-l\ almost all populations tested for these mecha- 
nisms were differentiated to species, with the majority of 
populations identified as An. gambiae s.s. (M and S form) 
and the remainder as An. arabiensis (Additional file 5). 
Relatively few An. gambiae s.l. populations were assessed 
for metabolic mechanisms via biochemical assays; when 
conducted, the majority of metabolic testing was for 
mono-oxygenases and populations were often not identi- 
fied to species. For An. funestus^ testing for metabolic 
mechanisms and modified acetylcholinesterase enzyme 
activity was more prevalent prior to 2001, with only two 
An. funestus populations tested for kdr mutations since 
2010. 

The prevalence of L1014S in An. gambiae s.l. popula- 
tions appeared to decrease after 2000 and then remain at 
moderate levels (39-56%) (Figure 5, Additional file 5). 
Conversely, the prevalence of L1014F populations ap- 
peared to increase after 2001-2003 and thereafter remain 
at higher levels (68 - 74%). While testing for metabolic 
mechanisms in An. gambiae s.l. was limited, the preva- 
lence of elevated activity of mono-oxygenases, esterases 
and GSTs appeared to decline between the periods 2004- 
2006 and 2007-2009, with mono-oxygenases and este- 
rases continuing to decline thereafter (to 24 and 43%, 



respectively) whereas detection of elevated GSTs became 
more common (67%). Testing for modified acetylcho- 
linesterase was limited to few assessments via enzyme as- 
says until methods for detecting the Ace-1^ mutation were 
developed, after which the gene was detected in 17 to 67% 
of An. gambiae s.l. populations depending on the period. 
No trends in the prevalence of resistance mechanisms 
were evident for An. funestus given the low number of 
populations tested. 

Distribution 

The spatial distributions between 2001 and 2012 of resist- 
ance to the four WHO-recommended adulticide classes 
and of the different resistance mechanisms are presented in 
Figures 6, 7, 8 and 9. These show that pyrethroid and DDT 
resistant Anopheles populations are widespread throughout 
the Afrotropical region, with foci of carbamate and organo- 
phosphate resistance reports in western Africa, Uganda, 
Ethiopia/northern Sudan and southern Africa. In terms of 
temporal trends in the geographic distribution of resistance, 
since 2001-2003 pyrethroid and DDT resistance have been 
reported in an increasing proportion of countries tested 
although this levelled off in 2010-2012 with resistance 
identified in 93% of the countries (Additional file 6). For 
An. funestus, populations have been tested in relatively few 
countries, confounding examination of the geographic 
spread of resistance. However, some spread in pyrethroid 
resistance is evident as prior to 2001 resistance was de- 
tected in one of four countries tested (Mozambique) 
whereas in 2010-2012 resistance was detected in the four 
countries tested (DRC, Kenya, Malawi, Mozambique). 
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Figure 5 Proportion of populations of An. gambiae s.l. with resistance mechanisms detected via standard biochemical or molecular 
methods. Note that populations are accounted for more than once when tested for multiple resistance mechanisms. Italics indicate the number 
of populations tested. 
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Both kdr mutations (L1014S and L1014F) are widely 
distributed in An. gambiae s.l. throughout the Afrotropical 
region (Figures 8 and 9). Metabolic mechanisms are also 
widely distributed although few populations from central 
Africa have been tested Ace-1^ has thus far been detected 
only in Anopheles populations from coastal areas of 
western Africa. Testing for resistance mechanisms in 
An.funestus s.s. has been sparse (Additional file 7). 

Tailored maps featuring these data and the additional 
PMI [34] and ANVR [33] datasets can be created via the 
IR Mapper online application (Figure 2). 

Discussion 

The number of scientific publications examining insecti- 
cide susceptibility and resistance mechanisms in Anopheles 
is growing rapidly. However, consolidated reviews of the 
status of resistance are few, and these can quickly become 
outdated due to the dynamic nature of insecticide resist- 
ance and the emergence of data for new or re-tested loca- 
lities, insecticides or vector species. Main resistance status 
reviews were published for Africa by the ANVR in 2005 
and 2012 [33,43], Coleman et al. in 2006 [9], Coetzee & 
Koekemoer in 2013 for An. funestus [44] and for pyre- 
throids by Ranson et al. in 2011 [10], with country profiles 
published by PMI [34] and an overview of recent (and 
often unpublished) data from malaria-endemic countries 
included in GPIRM [11]. These noted increasing resis- 
tance, but did not examine the full published historical 
data set for temporal and spatial trends. 

The database displayed via the online geospatial appli- 
cation IR Mapper consists of data extracted from public 



documents, with monthly additions capturing informa- 
tion released in new publications. Review of a subset of 
the data from evaluations conducted on Anopheles from 
the Afrotropical region according to standard proce- 
dures [26] indicated that while there are data available 
for 28 African countries, testing is often concentrated in 
particular sentinel or study sites where appropriately 
skilled entomologists operate. These may not represent 
the most epidemiologically significant zones but rather 
may be favoured due to their proximity to research faci- 
lities. Due to this heterogeneous - and for many areas, 
sparse - distribution of resistance data, and also owing 
to the highly focalised and dynamic nature of insecticide 
resistance, producing region-wide estimates of resistance 
status may have limited utility for programmatic pur- 
poses. The exception may be if clear drivers of resis- 
tance, such as agricultural usage of insecticides, can be 
characterised and extrapolated for specific regions. 

The high number of insecticide susceptibility reports for 
An. gambiae s.l. versus individual species of this complex is 
an indication of how seldom differentiation of sibling spe- 
cies is conducted post-bioassay. This is similarly the case 
for metabolic resistance mechanism assessments via bio- 
chemical assays. The paucity of data for individual mem- 
bers of the An. gambiae species complex and the An. 
funestus group limited comprehensive evaluations for 
these species and therefore this review focussed on wider 
species groupings. There is a critical need to identify spec- 
imens from susceptibility tests or mechanisms assays to 
species level, since insecticide susceptibility is clearly 
species-dependent. Further data are also required for 
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Figure 6 Distribution of insecticide resistance in An. gambiae s.l. collected between 2001 and 2012. Where there were multiple 
collections, species or insecticides tested, the lowest susceptibility category is displayed. Shading indicates the predicted distribution of the 
species complex [42]. 



An. funestus, which is a widespread and important vec- 
tor especially in central and southern African countries 
but has been under-tested likely due to the relative diffi- 
culty of rearing adults from the progeny of wild blood- 
fed females or larval field collections. The lack of data 
for certain regions or species may also be due to a limited 
local capacity to conduct simple bioassays as well as pro- 
pensity not to report data where widespread susceptibility 
is detected. 

Resistance mechanisms testing for An. gambiae has 
largely focussed on target site mutations {kdr and insensitive 



acetylcholinesterase G119S [Ace-1 ]) with relatively few 
resistant populations assessed for metabolic mechanisms. 
Evaluations of resistance mechanisms in An. funestus have 
been very sparse and mostly assessed metabolic me- 
chanisms and acetylcholinesterase enzyme activity. Urgent 
attention should be afforded to testing for resistance 
mechanisms particularly in An. funestus populations given 
the emerging evidence on the relative importance of this 
species as a primary vector in an increasing number of 
countries in the Afrotropical region [44,45]. In particular, 
the focus should be on metabolic mechanisms, as it is 
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Figure 7 Distribution of insecticide resistance in IKn. funestus s.l. collected between 2001 and 2012. Where there were multiple collections, 
species or insecticides tested, the lowest susceptibility category is displayed. Shading indicates the predicted distribution of the species 
group [42]. 



generally accepted that metabolic-based resistance is likely 
to have more severe implications than target-site resis- 
tance [11]. Such information is crucial to guide evidence- 
based insecticide resistance management and deployment 
of the most appropriate tools. 

The increasing prevalence and distribution of pyreth- 
roid and DDT resistance in An. gambiae si from the 
Afro tropical region is alarming. The GPIRM [11] was 
presented to deal with this, through the implementation 
of insecticide resistance management strategies aimed at 
preserving the efficacy of current tools. There is also a 



critical need to determine the impact that phenotypic 
resistance may be having on the efficacy of control tools 
across a range of eco-epidemiological settings, as well as 
factors contributing to such resistance and associations 
with the various resistance mechanisms. 

Given our ongoing reliance on insecticidal interven- 
tions, it is likely that the necessity and demand for sur- 
veillance data will increase especially as concerns mount 
on the impact resistance may have on malaria vector 
control. In order to support informed deployment of in- 
secticidal interventions, it is clear that resistance data 
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Figure 8 Distribution of molecular / biochemical resistance mechanisms in An. gambiae s.l. collected between 2001 and 2012. For sites 
for which multiple collections were tested, 'detected' is shown in preference to 'not detected'. Shading indicates the predicted distribution of the 
species complex [42]. 



must be collected, collated and fed back into the decision- 
making process in order to be optimally responsive to the 
local situation [26] . IR Mapper was primarily developed to 
collate essential data on insecticide resistance and make it 
readily available and easily digestible for malaria vector 



control programme managers and specialists. The tool 
allows rapid assessment of: (i) the geographic extent and 
frequency of resistance monitoring for specific vectors, 
insecticide classes, individual insecticides or mechanisms 
in a given region, and (ii) trends in resistance status and 
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Figure 9 Distribution of molecular / biochemical resistance mechanisms in An. funestus s.l. collected between 2001 and 2012. For sites 
for which multiple collections were tested, 'detected' is shown in preference to 'not detected'. Shading indicates the predicted distribution of the 
species group [42]. 
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detected mechanisms over specified periods. Such in- 
formation is used to identify major gaps in insecticide 
susceptibility and resistance mechanisms monitoring, as 
well as to ascertain the most up-to-date status based on 
published information. With the addition of users' own 
data, unpublished data sets can also be included in such 
evaluations. 

In terms of guiding interventions, while IR Mapper can- 
not directly inform product deployment, which depends 
on many other variables including price, availability, and 
suitability for intended end users, it does have a clear role 
in informing the selection of appropriate classes or insec- 
ticides for use in IRS as part of an evidence-driven insecti- 
cide resistance management strategy. The ability of IR 
Mapper to zoom in to targeted districts within a specific 
country provides the means for making large amounts of 
data far more easily comprehensible as well as highlighting 
regions where data are lacking. These visual maps are in- 
valuable for monitoring trends in resistance development 
over time and planning insecticide-based control interven- 
tions accordingly. For instance, high resistance to pyre- 
throids (e.g. mortality in WHO susceptibility test of <50%) 
and a rapid increase in the frequency of kdr mutations in 
an area earmarked for IRS will suggest that an immediate 
switch to an acetylcholinesterase inhibitor such as a 
carbamate or organophosphate is justified as per 
GPIRM technical recommendations [11]. Alternatively, 
if susceptibility remains and kdr mutations are absent 
or frequency is stable, it may indicate that pyrethroids 
can continue to be used if applied pragmatically in con- 
sideration of the principals of pre-emptive resistance 
management. IR Mapper may also be used to advocate 
for increased monitoring in regions where specific 
issues of high resistance have been identified and also 
to determine if there are any changes over time in either 
insecticide susceptibility or resistance gene frequency 
following specific interventions. 

Data contained on the IR Mapper platform should be 
interpreted with consideration of the limitations of the 
methods applied to derive the displayed data. WHO in- 
secticide susceptibility and resistance mechanisms data 
cannot be used to draw conclusions about possible oper- 
ational control failures of insecticidal vector control tools 
such as LLINs or IRS. For this, comprehensive field evalu- 
ations correlated with malaria case data would need to be 
conducted. Moreover, since IR Mapper does not allow for 
filtering of data on the basis of assay parameters such as 
insecticidal dosage on impregnated test papers, or the 
number of mosquitoes tested, users should be cautious in 
interpreting the data en mass. Relevant information is pre- 
sented for each individual data point as is an online link to 
the original data source, and the onus is on the user to as- 
certain the relative importance of available data for guid- 
ing evidence-based surveillance and control strategies. 



Although IR Mapper is currently limited to data from 
public scientific publications and a few other credible data 
sources (e.g. PMI, ANVR), incorporation of unpublished 
data sets is feasible if a robust verification and authorisa- 
tion process in established. This would require a sustained 
input of labour considering the increase in volume of in- 
secticide resistance data becoming available. Results from 
microarray studies or for newly-identified resistance 
mechanisms (such as new sodium channel mutations) can 
easily be incorporated on the existing configuration, since 
the user interface automatically reflects the parameters 
included in the database as long as the outputs of 
detected/not detected are defined. There is also the possi- 
bility to extend the platform for display of CDC bottle bio- 
assay data [46], which may be important considering this 
is now being applied as the primary insecticide suscepti- 
bility surveillance method in some countries. 

Conclusions 

IR Mapper is a simple tool for investigating temporal and 
spatial trends in Anopheles insecticide susceptibility and re- 
sistance mechanisms. National malaria control programs 
can use the platform to gain an overview of the spatial dis- 
tribution and extent of available data to inform surveillance 
strategies. Outputs can be used to optimise malaria vector 
control via pragmatic use of insecticidal interventions and 
design of insecticide resistance management strategies. It is 
critical that such evidence-based approaches are prioritised 
given the apparent trend of increasing pyrethroid and 
DDT resistance in Anopheles in the Afrotropical region. 

Additional files 



Additional file 1: Number of populations for which resistance was 
confirmed and not confirmed using standard WHO insecticide 
susceptibility tests on adult mosquitoes. A) An. gambiae s.l.; B) An. 
gambiae s.s.; C) An. arabiensis; and D) An. funestus s.l. by year of mosquito 
collection and insecticide class. Insecticide susceptibility testing was most 
common for /\n. gannbiae s.l. populations. Pyrethroids and 
organochlorines were the most frequently tested insecticide classes. 

Additional file 2: Number of populations for which resistance 
mechanisms were detected and not detected. A) An. gambiae s.l., 
B) An. gannbiae s.s., C) An. arabiensis and D) An. funestus s.l. by year of 
mosquito collection and mechanism class. An. gannbiae s.l. was 
commonly tested for mechanisms while An. funestus s.l. were not. kdr 
mutations were the most frequently tested resistance mechanisms in An. 
gannbiae s.l. with metabolic mechanisms seldom tested. 

Additional file 3: Number of Anopheles populations for which 
insecticide susceptibility and resistance mechanisms tests were 
conducted between 2001 and 2012 for top 20 malaria burden 
countries. Few or no Anopheles populations were tested for insecticide 
susceptibility or resistance mechanisms in some of the countries with the 
highest malaria burden. 

Additional file 4: Summary of data from WHO insecticide 
susceptibility tests conducted with Anopheles spp. populations 
collected in Africa between 1963 and 2012. Prevalence of confirmed 
pyrethroid and organochlorine resistance in An. gannbiae s.l. increased 
over time. 
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Additional file 5: Summary of data from resistance mechanisms 
tests conducted with Anopheles spp. populations collected in Africa 
between 1987 and 2012. Prevalence of L1014S remained at moderate 
levels after 2001-2003 whereas prevalence of L1014F was at 
comparatively higher levels. 

Additional file 6: Number of countries in Africa for which resistance 
was confirmed in at least one population of Anopheles spp. via 
WHO insecticide susceptibility tests. Since 2001-2003, pyrethroid and 
DDT resistance have been reported in an increasing proportion of 
countries tested although this levelled off in 2010-2012. 

Additional file 7: Number of countries in Africa for which resistance 
mechanisms were detected in at least one population of Anopheles 
spp. Testing for resistance mechanisms in An. funestus s.s. has been 
sparse. 
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